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Abstract 

The decays of the Z 6 (10610) and the Z fc (10650) to T(nS>, h b (mP)n and Xbj(mP)l 
(n = 1,2,3, m = 1,2 and J = 0,1,2) are investigated within a nonrelativistic effective 
field theory. It is argued that, while the decays to Y(nS)ir suffer from potentially large 
higher order corrections, the P-wave transitions of the Zj, states are dominated by a single 
one loop diagram and therefore offer the best possibility to confirm the nature of the 
states as molecular states and to further study their properties. We give nontrivial and 
parameter-free predictions for the ratios of various partial widths of the Z& and Z' h into 
final states with hb(mP)TT and Xbji^P)! ■ While such relations appear naturally in the 
molecular picture for the mentioned transitions, they are not expected to hold for any 
other scenario. In addition, the branching fractions for the neutral ^-states to XfcJT are 
predicted to be of order 10 -4 -10 -3 . This provides a fine test of the molecular nature in 
future high-luminosity experiments. 
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1 Introduction 



Recently the Belle Collaboration found two narrow structures, namely the ^(10610) = Z b 
and the Z£(10650) = Z^, in the T(n5)vr ± (n = 1,2,3) and k b (mP)ir L (m = 1,2) invariant 
masses of the T(5S) — > T(nS)7r + 7r _ and Y(5S) — > hb(mP)ir + Tr~ decays processes [I]. In 
their latest data, also the open-bottom channels B*B and B*B* are seen [3J. The fact that 
they lie in the bottomonium mass region and they are charged means that they cannot be the 
conventional bottomonium mesons, and their isospins are 1. The observation of the neutral 
state with a mass consistent with that of the .^(10610) [3] presents a further confirmation 
that they are members of isotriplets. If they exist as observed by the Belle Collaboration, 
they must be exotic states with a pair of hidden bb and valence light quarks. Considering 
parity and charge parity, their quantum numbers should be I G (J P ) = 1 + (1 + ), and the charge 
parity of the neutral state is negative. 

Assuming that the total width of Z^' is saturated by the seven channels already observed 
experimentally, i.e. T(nS)vr (n = 1,2,3), h b (mP)7r (m = 1,2), BB* + B*B and B*B*, Belle 
gives the branching ratio of each channel in Y(5S) three-body decays [2j. The proximity 
of the states to the B^B* thresholds leads to the suggestion that they could be hadronic 
molecules of the corresponding states [U El El El El El [101 EEU EE2], to be distinguished from 
the compact bqbq tetraquarks [HJ Q31 [2] . By hadronic molecules, we mean states composed 
of hadrons — they can be bound states (poles on the physical sheet with respect to the 
bottom- meson channel with a mass smaller than the threshold value), resonances or virtual 
states (both on the second sheet with respect to the relevant bottom-meson channel — the 
former above, the latter below the threshold). Based on a nonrelativistic effective field theory 
(NREFT) HSldB], the authors of Ref. [5] show that the /i 6 (lP, 2P)7r ± data can be described 
within the bound state scenario. Therein, the Z^B^B* coupling constants are related to the 
binding energies using a model-independent relation for S-wave shallow bound states |171I18], 

However, since now data on the decays of z9 to open bottom channels are available, we can 
chose a more general ansatz and take these couplings from data directly — in this way our 
results are valid for bound states, resonances as well as virtual states. Since we start from 
the assumption that the Z b states are purely molecular states, their decays into the T(nS)ir 
and hb(mP)ir can only happen via BB* + B*B (for simplicity, BB* will be used to represent 
BB* + B*B in the following) and B*B* loops. Since both Z^s are located very close to 
the corresponding open-bottom threshold, the system can in principle be examined by the 
NREFT approach. 

As shown in Refs. [151 H6j . a systematic power counting can be established for the ir or r\ 
emissions between charmonium states. Because the S'-wave and P-wave heavy quarkonia cou- 
ple to the open-flavor heavy meson and anti-meson in a P-wave and an S'-wave, respectively, 
the transitions studied in Ref. [16] can be classified into three groups, namely transitions be- 
tween the S-wave heavy quarkonium states, P-wave states, and between the P- and S-wave 
states. The decay amplitudes of different groups have their own nonrelativistic velocity count- 
ing. Because the Z b states have positive parity, they couple to the bottom and anti-bottom 
mesons in an S-wave. Thus, the transitions Z b — > T(nS)7r are analogous to those between 
P- and S-wave quarkonia, while the Z\> — > hb(mP)ir processes are similar to those between 
two P-wave quarkonia. The main difference is that the normal heavy quarkonium transitions 
with the emission of a pion studied in Ref. [16] break isospin symmetry while the Z\, decays 
do not. Thus, by studying the (Z' b ) — > T(nS)ir and /i&(mP)7r, we can examine the power 
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counting rules established in Ref. [16], and also better understand the properties of these two 
exotic states. 

Further insight can be gained by studying radiative decays of the neutral Z& states into 
bottomonia. Because the Xbj (J = 0, 1, 2) states are the spin partners of the hb of the same 
principal quantum number, they couple to the bottom mesons with the same coupling constant 
in the heavy quark limit. Thus, since molecular states can decay via bottom meson loops only, 
the radiative decays Z b — > Xbj(nP)"f are related to the pionic decays zjf — > h^nP)^. 
Note that these transitions would be unrelated if the Z b states were of tetraquark nature. 
Similar considerations were made in Ref. [19] for hindered Ml transitions of the P-wave 
charmonia. These decay channels have not been observed so far, but they could be potentially 
important in confirming the molecular nature of the Zb states, and thus are worthwhile to 
study experimentally. 

In this paper, we will assume that the Z^ are dynamically generated from the B^B* 
interactions, i.e. hadronic molecules of the B^B*. We will try to identify the quantities 
which are sensitive to such a scenario. Section [2] contributes to the power counting in the 
NREFT framework of the relevant decays, which are the decays of the Zb states into the Tit, 
hb 71 ', XbJl — in particular we show that not all decays are accessible to the formalism. The 
numerical results are given in Sec. El In SecHJ we compare our results to previous calculations 
and make some comments. A brief summary is presented in Sec. [5] The loop function used 
in the calculations and the decay amplitudes in the NREFT are collected in Appendix lAl As 
a cross check, we also calculate the same quantities using a Lorentz covariant formalism, and 
the formulas are summarized in Appendix [Bl 

2 Power Counting 

In Ref. [15j . a NREFT method was introduced to study the meson loop effects in the heavy 
quarkonium transitions. The power counting scheme was analyzed in detail in Ref. [16]. The 
key quantities here are the velocities of the intermediate mesons. In this section, we briefly 
review the ideas of Refs. [151 116] together with an improved discussion for the higher loop 
diagrams. 

In general, the heavy meson velocities relevant for the decay of some particle X may be 
estimated as vx ~ y\Mx — 2mB\/mB, where the absolute value indicates that the formula 
can be used for both bound systems as well as resonances. The analogous formula holds when 
the two heavy mesons merge to a quarkonium in the final state. According to the rules of 
a nonrelativistic effective field theory [20] (for a review, see e.g. [21]), the momentum and 
non-relativistic energy count as vx and v x , respectively. For the integral measure one finds 
v x /(4tt) 2 . The heavy meson propagator counts as ^jv\. The leading order 5- wave vertices do 
not have any velocity dependence, while the case for the P-wave vertices is more complicated: 
it scales either as vx when the momentum due to P-wave coupling contracts with another 
internal momentum, or as the external momentum q when q is contracted. 

We start with the radiative transitions as shown in the upper row of Fig. [TJ If the Zb states 
are molecular states, their spin wave functions contain both s bb = and 1 components @], 
where s bb is the total spin of the bb component. Thus, the radiative decays of the Zb states 
into the spin-triplet Xbj can occur without heavy quark spin flip and survive in the heavy 
quark limit. This is different from the Ml transitions between two P-wave heavy quarkonia 
which have been analyzed in Ref. [19] . Both the couplings of Zb and Xbj to a pair of heavy 
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Figure 1: Schematic one- and two-loop diagrams of the transitions Zb —> Tir, h^ir and XbJl- 



mesons are in an S— wave, and the photon coupling to the bottom mesons is proportional to 
the photon energy EL. For the diagram of Fig. [1] (a) the amplitude therefore scales as 
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(1) 



where the velocity that appears is v = [yz + v x )/2 ~ u x /2 [22], since ~ u^/ ~ 0.02, if the 
central values of the measured Z^} masses, 10607.2 MeV and 10652.2 MeV, are used, while 



v x ranges from 0.12 for the Xbj(3P) to 0.26 for the Xbj(2P) to 0.37 for the XbjO-P)- Here we 
used the mass of the x&j(3P), 10.53 GeV, as reported by the ATLAS Collaboration [23]. In 
the following, we will count v as 0(v x ). The scaling ensures that the amplitude gets larger 
when the bottomonium in the final state is closer to the open-bottom threshold. Thus, we 
expect for the absolute value of the decay amplitude from this diagram 
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: : = 1 : 1.4 : 3.1, 



(2) 



if the Xbj{ n P)BB coupling constants take the same value. Diagram (a) can be controlled 
easily in theory. Thus, clear predictions can be made whenever diagram (a) dominates. In the 
following we will identify such dominant decays based on the power counting for the NREFT. 

As for Fig. [T] (b), the coupling XbBB^f cannot be deduced by gauging the coupling of Xbj 
to a BB-meson pair. Thus, it has to be gauge invariant by itself and proportional to the 
electromagnetic field strength tensor F^ v . This gives a factor of photon energy and the 
amplitude of Fig. [T] (b) scales as 



(4vr) 2 (wfvA 



E^V Z 



(3) 



where we have assumed that the corresponding coupling is of natural size. Thus, diagram (b) 
is suppressed compared to diagram (a) at least by a factor of vzv x < 0.01 for the decay to 
Xfej(l-P) and even smaller for the excited states. 

The situation is more complicated for the graph displayed in Fig. [T] (c). Here we have a 
two-loop diagram, so that the velocities running in different loops are significantly different 
- the one in the loop connected to the Zb is vz, and the other is v x . It is important to count 
them separately since vz *C v x due to the very close proximity of the Zb to the threshold^. The 

lr The concept applied here is analogous to the scheme by now well established for the effective field theory 
for reactions of the type iVJV — > NNtt — see Ref. [U] for a review. 
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internal pion momentum scales as the larger loop momentum, and thus the pion propagator 
should be ~ l/(m B v 2 ). This leaves us with 

4 1 v x 1 1 E i9 9 4 vz E ig 2 m 2 B 

B „. / / i_\2a2 ' \ ' 



(47r) 2 (4) 2 (ATT) 2 («2)2 m |t;2 ^ P w * % ^A 2 , 

where P^ is the pion decay constant in the chiral limit, the factor m B has been introduced 
to give the same dimension as the estimate for the first two diagrams, and the hadronic scale 
was introduced via the identification A x = 4:irF n . Thus the two- loop diagram is suppressed 
compared to the leading one, Eq. ©, by a factor vzg 2 m B / A 2 , ~ 0.1, where we used for the 
coupling B* — > Btt the value g = 0.5 (a recent lattice calculation gives 0.449 ± 0.051 [25]). 
and A x ~ 1 GeV. It can easily be seen that Fig. [1] (d) gives the same contribution which also 
reflects the fact that they are both required at the same order to ensure gauge invariance. 
Thus, from our power counting it follows that the loop diagrams of Fig. [JJ (b)-(d) provide 
a correction of at most 10%. We will therefore only calculate diagram (a) explicitly and 
introduce a 10% uncertainty for the amplitudes which corresponds to 20% for the branching 
ratios. Higher loop contributions are to be discussed later. 

Next we consider the hadronic transition Zb — > h^n. This decay has already been studied 
in Ref. [5] in the same formalism. Again, since the has even parity, its coupling to the 
bottom mesons is in an S-wave. In addition, the final state must be in a P-wave to conserve 
parity, such that the amplitude must be linear in the momentum of the outgoing pion, q. We 
therefore find for the one-loop contribution of Fig. [T] (e) 

v 5 1 gq qF n 

9^ttt , (5) 



(4tt) 2 (u 2 ) 3 F„ vA\ 
while Fig. [1] (f) gives 



(6) 



(4^)2(4)2^ A| • 

Notice that the pion has to be emitted after the loop if the Z^ is a pure hadronic molecule, 
so that the velocity in the counting should be vz instead of Vh- This is suppressed by VhVz/g 
which leads to a correction of the order of 2% noticing that ~ v x . 
The two-loop diagram Fig. [T] (g) contributes as 



1 v x 1 1 gqE* 3 v z gF v q z m B 



m i B = --^r^, (7) 



(4vr)2(4)2( 4 7r)2(4)2 m |4^F2 " V h A\ 

where we have used that the energy from the ttB —> ttB vertex can be identified with the 
energy of the outgoing pion [26] . E n ~ q. The B*Bir vertex contributes a factor of the 
external momentum q since the Z\, —> h^ir is a P-wave decay, and this is the only P-wave 
vertex. Therefore, this diagram is suppressed compared to the leading loop, Eq. ([5]), by 
a factor vzrriBq/ A x which is smaller than 10%. Thus, also for the transitions Z^h^-K we 
may only calculate the leading one- loop diagrams, Fig. [T] (e), and assign an uncertainty of 
10% to the rates which gives an uncertainty of 20% for the branching ratios. Higher loop 
contributions are to be discussed later. 

Finally, we consider at the decay channel Z^ — > Ttt. Here the final state is in an S-wave, 
but the coupling of the T to BB is in a P-wave. For diagram (e) the momentum due to this 
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Figure 2: Two three-loop diagrams contributing to the decays of the Zj, into a heavy quarko- 
nium and a pion or photon. 

coupling has to scale as the external pion momentum. Together with the pionic coupling that 
is also linear in the pion momentum, the amplitude is thus proportional to q 2 . The one-loop 
diagram for Z b — > Ttt via Fig. [T] (e), is therefore estimated as 

v 5 q gq 



(4tt) 2 (v 2 ) 3 Ft, y vK\ 
The diagram Fig. Q] (f) on the other hand gives 



1 E n vzE^F^mB 

-m B ~ tt; , (9) 



(4vr)2(4)2^ ° 1\\ 

where wig is introduced in order to get the same dimension as Eq. ([8]). Compared to the 
one- loop diagram (e) this is a relative suppression of order vrvzniB/q which is less than 10% 
for all the T states, where the value of t>x is about 0.46, 0.33 and 0.22 for the IS", 25 and 
3S states, respectively. The two- loop contribution with the exchange of a pion, Fig. CD (g), is 
estimated as 

v% 1 v\ 1 1 v\g E w 5 vzVygE^ 3 

T m B 71 m B ■ ( 10 ) 



(4tt) 2 (4) 2 (4vr) 2 (4)2 m 2 B v\ F n F 2 B k\ 

Thus, the strength of two-loop diagram relative to leading one- loop diagram is estimated for 
the Ttt as v zv\m B / '(A^q) . Numerically, this corresponds to a factor of around 0.6 for the 
T(1S,2S)tt and 0.7 for the T(3S')7r amplitudes. As a consequence, the branching ratios for 
these transitions can only be calculated with rather large uncertainties up to 100%. 

The heavy meson velocities relevant for the mentioned transitions range from 0.02 to 0.5 
- in momenta this is a range from 0.1 to 2.5 GeV. While pion contributions are expected to 
be suppressed significantly and can be controlled within chiral perturbation theory for pion 
momenta of up to 500 MeV and smaller, higher pion loop contributions might get significant 
for momenta beyond 1 GeV. We will now study those higher pion loops within the power 
counting scheme outlined above. We will start with the three-loop diagrams, as shown in 
Fig. [2] considering first diagram (a). The results can be easily generalized to higher loops 
as shown below. Compared to the two-loop diagrams in Fig. [TJ there are one more pion 
propagator and two more bottom meson propagators in the transition and none of them 
is connected to the external heavy quarkonia. In addition, there is no two-bottom-meson 
unitary cut present. As a consequence, we have to use a relativistic power counting — c.f. 
Ref. |20] , Then pion momentum and energy are of order iriBV b i, with v b i the velocity of the 
5-meson connected to the 66-meson in the final state. Both the energies and momenta of the 
additional bottom mesons are now of the same order, such that the bottom-meson propagator 
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Table 1: Left: Preliminary measurements of the branching ratios for zj~ from the Belle 
Collaboration i2j. Right: Masses of the various particles used here [27] [28] . 



Branching ratio (%) 


Zj, (10610) 


Z' b (10650) 


T(15)vr+ 


0.32 ± 0.09 


0.24 ± 0.07 


T(25)tt+ 


4.38 ± 1.21 


2.40 ± 0.63 


T(35)vr+ 


2.15 ±0.56 


1.64 ± 0.40 


h b (lP)ir+ 


2.81 ± 1.10 


7.43 ± 2.70 


h b (2P)w+ 


4.34 ± 2.07 


14.82 ±6.22 


B+B*° + B°B*+ 


86.0 ±3.6 




B*+B*° 




73.4 ± 7.0 





Mass [GeV] 




Mass [GeV] 


T(1S) 


9.460 


Xbo(lP) 


9.859 


T(2S) 


10.023 


Xbi(lP) 


9.893 


T(3S) 


10.355 


Xw(l-P) 


9.912 


h b (lP) 


9.899 


XM)(2P) 


10.233 


h b (2P) 


10.260 


Xbl(?P) 


10.255 


B 


5.279 


Xw(2P) 


10.269 


B* 


5.325 


7T 


0.138 



is counted as l/v bb . The integral measure reads v^/ (An) 2 . Therefore, the additional factor 
as compared to the two-loop diagrams is 



v 



4 -i / \2 



1 (9 v bb) _2 _ ( 9m B v h i 



If v bb ~ 0.4, then the three-loop diagram (a) is of similar size as that of the two-loop diagrams. 
This is the case for the processes with the hb(lP), x&j(1P), and Y(1S, 2S). For smaller values 
of v bb , it is suppressed. In diagram (b) there is only one more bottom meson propagator and 
the additional BttBtt vertex is in an S"-wave. We obtain 

Vtl 1 Vb \rnl=(^P*) 2 . (12) 



(4vr)^3_ F 2 « V A x y 

Four and higher loop diagrams that cannot be absorbed by using physical parameters may 
now be estimated by applying a proper number of factors of the kind of Eqs. (jlip and (|12p . 
It is easy to see that also additional topologies provide analogous factors. Since mg/Aj, ~ 5, 
higher loops get increasingly important, if v bb > 0.2. For v bb ~ 0.2, the three and more loops 
are of the same order as the two-loop diagrams, which is the case for the 2P states, and thus 
suppressed in comparison with the one-loop contribution. For the 3P bottomonia in the final 
state, the value of v bb is even smaller, and the multiple loops are even suppressed as compared 
to the two-loop contribution. 

To summarize the findings of the power counting analysis, we conclude that the calculation 
of one-loop triangle diagrams as depicted in Fig. Q] (a) and (e) is a good approximation, with 
a controlled uncertainty, to the transitions of the Z b into the Xfcj(2-P, 3P)7 and h b (2P)-n - 
for the h b (3P)ir the phase space is too limited. But similar calculations are not applicable 
to the decays into the IP states as well as the T(nS)ir. For the decays Z b — > T(nS)n, the 
contribution from the two-loop diagrams is not suppressed or even enhanced compared to 
the one-loop contribution. In light of this discussion, it becomes clear why the pattern of 
branching fractions for these channels (see Tab. CE} £[T(2S>] > B[T(3S)tt] > B[T(13>] 
cannot be reproduced by calculating the three-point diagrams in the NREFT formalism, 
which always favors the IS" to the 2S transition and the 2S compared to the 3S transition 
due to the factor q 2 in Eq. (|SJ). 
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3 Results 



In this section, we investigate quantitatively the decays of the Z b states through heavy meson 
loops. The coupling of these states to the heavy meson fields H a = V a • a + P a and H a = 
—V a ■ a + P a with V a (Va) and P a (P a ) annihilating the vector and pseudoscalar (anti-)heavy 
mesons, respectively, is given by the Lagrangian 



£ z = i-{Zl l a H a a l H b ) + R.c. 
where the Z b states are given by a 2 x 2 matrix 




Zla z l —Z 0i 



z +i 



ba 



We incorporate the experimental observation that the ZjJ couples only to B^B* via the 
Lagrangian, 

C zz , = z > £ vky\i Z j V \k + z \yV Z ip\ _ p\z i V^ i \ + H.C. 
Fitting to the experimental data we find 

z = (0.79 ± 0.05) GeV" 1/2 , z' = (0.62 ± 0.07) GeV" 1/2 , (13) 

and especially 

J = 1.27 ±0.16 (14) 

which deviates from unity by 2a. This deviation indicates a significant amount of spin sym- 
metry violation, which, however, is not unnatural for very-near-threshold states where small 
differences in masses may imply huge differences in binding energies resulting in significantly 
different effective couplings, as already discussed in Ref. [5]. The nonrelativistic Lagrangian 
for the Xbj coupling to a pair of heavy mesons can be found in Ref. [29] , and the one for the 
magnetic coupling of heavy mesons in Ref. [30] . 

With the amplitudes given in Appendix [Al it is straightforward to calculate the de- 
cay widths of the zj: — > hb(mP)TT, which are proportional to gf, where g\ is the P-wave 
bottomonium-bottom meson coupling constant. In the ratio defined as 

_ T[Z' b ^ h b (mP)n] 
U - T[Z b h b (mP)nY [Lb) 

gi is cancelled out. With the meson masses listed in Tab. Q3 we obtain 



£i = 1-21 



0.75, & = (1-53 ± 0.43) 



0.95 ±0.36. (16) 



where the first error in the second term is the theoretical uncertainty due to neglecting higher 
order contributions (see the discussion in the previous section), and the second one also 
includes the uncertainty of z' / z added in quadrature. Due to the theoretically uncontrollable 
higher order contributions for the decays into the IP states, no uncertainty is given for £i. 
The predictions are consistent with their experimental counterparts (see Tab. [1]) 

£f x P = 1.65 ± 0.96, £2 xp = 2.13 ± 1.44. (17) 
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Table 2: The ratios £ of different decay modes in both NREFT and a relativistic framework 
as compared with the experimental data. The NREFT values quoted without uncertainties 
may be understood as order-of-magnitude estimates. 





NREFT 


Rel. 


Exp. 


T(1S)tt 


0.7 


0.7 


0.47 ±0.22 


T(25)vr 


0.9 


0.8 


0.34 ±0.15 


T(35)vr 


2±2 


1.6 


0.48 ± 0.20 


h b (lP)n 


0.8 


0.7 


1.65 ±0.96 


h b (2P)7T 


1.0 ±0.4 


0.9 


2.13 ± 1.44 



Here, new measurements with significantly reduced uncertainties would be very desirable. A 
collection of ratios for the decays of the Z b states to h^mP)^ and Y(nS)n, respectively, is 
presented in Tab.[2j The uncertainties, whenever they are under control theoretically, are also 
included. The significant deviations for the T(nS)7r results from the experimental numbers 
appear natural, given that for those transitions higher loops were argued to be at least as 
important as the one-loop diagram included here, as outlined in detail in the previous section. 
As a cross check of our nonrelativistic treatment, we also calculated the same quantities using 
a Lorentz covariant formalism with relativistic propagators for all the intermediate mesons 
(the formalism is summarized in App. [B]). The results are denoted as 'Rel.' in Tab. [2j They 
should be compared with the central values of the NREFT results. The difference reflects 
relativistic corrections to the nonrelativistic treatment of the bottom-meson propagators. In 
our case, the difference between relativistic and NREFT calculations here and in the following 
never exceeds 15% for the rates and thus is well below the uncertainty due to higher loops. 

It is important to ask to what extent the above predictions can be used to probe the 
nature of the Z b states. The Z b and Z' b are away from the BB* and B*B* thresholds by 
similar distances, mz b — tub — itib* — m z' b ~ 2ms*. Additionally, due to the heavy quark 
spin symmetry, one may expect that the loops contribute similarly to the decays of these two 
states into the same final state. This is indeed the case. We find that the ratios aside of 
\z' / z\ 2 are basically the phase space ratios, which are 

m -» Miwi 3 _ L20i m -> wpmi' = 1 53 (18) 



\q(Z b -> /i,,(1P)i)| 3 ' \q(Z h -> M2-PM 

This implies that the ratios for the decays into the hf,(mP)ir are determined by the ratio of 
the partial decay widths for the open-bottom decay modes, 

PS' m T(Z'+ -> B*+B*°) 
?m PS m T(Z b + — >• B+B° ± B°B* + ) ' 1 ; 

where PSm is the phase space for the decays Z^ — > hb{mP)iT. Such a relation cannot be 
obtained were the Z b states of tetraquark structure because the decay of a bqbq tetraquark 
into the h^-rr knows nothing about the decay into the open-bottom channels. Here, if we 
impose spin symmetry for z and z', one would get for £ m simply the ratio of phase spaces. In 
reality the second factor which is the square of the ratio of effective couplings, c.f. Eq. (|14p . 
deviates from unity due to spin symmetry violations enhanced by the proximity of the B^B* 
thresholds as discussed above. 
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Table 3: The ratio f2 and the corresponding branching fractions for all possible radiative 
decays. Uncertainties are given, whenever they can be controlled theoretically (see text). 
The values quoted without uncertainties may be understood as order of magnitude estimates. 







z b 




K 




n 


Branching Fraction 


n 


Branching Fraction 




5 x 10" 3 


1 x 10" 4 


4 x 10" 3 


3 x 10" 4 


X6i(lP)7 


i x icr 2 


3 x 10~ 4 


1 x 10~ 2 


8 x 10" 4 


X62(1P)7 


2 x 10~ 2 


5 x 10~ 4 


2 x 10~ 2 


1 x 10~ 3 


X6o(2P)7 


(6.3 ± 1.8) x 10~ 3 


(2.7 ± 1.5) x 10~ 4 


(4.2 ± 1.2) x 10~ 3 


(6.2 ±3.2) x 10~ 4 


X6i(2P)7 


(1.3 ±0.4) x icr 2 


(5.6 ±3.2) x 10~ 4 


(1.3 ±0.4) x 10~ 2 


(1.9 ± 1.0) x 10~ 3 


Xb2(2Ph 


(1.9 ±0.5) x 10" 2 


(8.3 ±4.5) x 10" 4 


(1.8 ±0.5) x 10" 2 


(2.7 ±1.3) x 10" 3 



Heavy quark spin symmetry allows one to gain more insight into the molecular structure. 
Because the x&j(mP) are the spin-multiplet partners of the hb(mP), the radiative decays of 
the neutral into Xbj{ m P)l can be related to the hadronic decays of the \ no matter 
whether they are neutral or charged, into the hb(mP)ir. It is therefore useful to define the 
following ratios 

n4' )0 XbAmPh ) 
T(ZP° -> h b (mP)n° 

With the coupling constants in the bottom meson-photon Lagrangian determined from else- 
where, see for instance Ref. [30] . such ratios can be predicted with no free parameters. At 
the hadronic level, the Z b radiative transitions can only be related to the hadronic ones if the 
Z b s are hadronic molecules so that the two different types of transitions involve the same 
set of coupling constants (modulo the bottom meson-pion/photon coupling which can be 
determined from other processes or lattice simulations). Thus, if the branching fractions of 
the radiative transitions are large enough to be detected, such a measurement would provide 
valuable information on the nature of the Zf, states. The results for the ratios are collected 
in Tab. [3j 

Using the branching ratio B(Z^ — > hi,(lP, 2P)ir Jr ), we find that the branching ratios 

of Z^jp — > Xbj(lP,2P)l are of order 10~ 4 ~ 10~ 3 . The largest branching fractions of the 
Xbj{ m P) are those into the jT(nS), and the T(nS) can be easily measured. Thus, the final 
states for measuring the Z® — > XbJl would be the same as those of the Z® — > T(nS)7r° be- 
cause the 7T° events are selected from photon pairs. This means that the detection efficiency 
and background of these two processes would be similar. In the preliminary experimental 
results [3], the Zf ) (W610)° event number collected in the T(1S', 2S')7r channels is of order 
0(100). Given that the luminosity of the future Super-KEKB could be two orders of magni- 
tude higher than KEKB, such transitions will hopefully be measured. Furthermore, one may 
also expect to measure these radiative transitions at the LHCb. Note that the experimen- 
tal confirmation of the ratios given in Tab. [3] would be a highly nontrivial evidence for the 
molecular nature of the Zf, states. 

Lacking knowledge of the Xbj{3P)BB coupling constant, the transitions into the 3P states 
cannot be predicted parameter-free. However, they can be used to check the pattern in Eq. ([2]) 
predicted by the power counting analysis. The decay widths of the Z b — > Xbj(n f iP)l are 
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Table 4: The ratios defined in Eq. ([2ip for all channels. Uncertainties are given, whenever 
they can be controlled theoretically. 



(J = ):(J = 1):(J = 2) 


Z h XbAlPh 1 = 2-5 
Z' b -> Xw(lP)7 1 = 2-9 
Z 6 ^Xbj(2P)7 1: (2.1 ±0.6) 
Z£ -> xw(2P)7 1 : (3.0 ± 0.9) 


3.7 
4.4 

(2.9 ±0.8) 
(4.2 ± 1.2) 



proportional to g\ mP . Taking the same value for g\ mP , the explicit evaluation of the triangle 
loops gives 1 : 1.8 : 4.4 for the ratios defined in Eq. ((2J) with J = 1. The values are close to 
the ones in Eq. ((2j), and thus confirm the 1/v scaling in Eq. ((TJ) of the amplitudes. 

As observed in Ref. |16j . the NREFT leading loop calculation preserves the heavy quark 
spin structure. Because the Z b contains both s bb = and 1 components, the leading con- 
tribution to its transitions into the normal bottomonia, which are eigenstates of s bb , comply 
with the spin symmetry. This conclusion should be true no matter what nature the Z b s have 
as long as the spin structure does not change. Thus, one expects that the branching fraction 
ratios of the decays of the same Z b into a spin-multiplet bottomonia plus a pion or photon, 
such as 

B(ZP -> Xboi^Ph) : B(Z 6 W -> XbiimPh) ■ B(Z® -> x 62 (mP) 7 ), (21) 

are insensitive to the structure of the Z^. This statement may be confirmed by observing 
that our results, as shown in Tab. [4] agree with the ratios 1 : 2.6 : 4.1 (for the IP states) and 
1 : 2.5 : 3.8 (for the 2P states) which are obtained solely based on heavy quark spin symmetry 
in Ref. |31j . One may expect a derivation from the spin symmetry results, which is due to the 
mass difference between the B and B* mesons, to be of order 0(2AQcr>/ m B) ~ 10%. The 
central values given in Tab. 0] deviate from the spin symmetry results by at most 20%, and 
they are fully consistent considering the uncertainties. 

4 Comparison with other works on decays 

Since their discovery in an impressive number of theoretical works the molecular nature of 
the Z b states was investigated. In this section we compare in some detail our approach to the 
calculations in Refs. [32] [33] [MIES] which deal with some of the decays considered in our paper. 
Common to most of these works is that, contrary to our approach, the second part of the one 
loop integral shown in Fig. Q] (e) is either approximated [32], [33] or calculated differently [M] . 
Especially, by this the analytic structure of the loop gets changed converting it to a topology 
of type (f) in that figure, since the second £?(*)_B(*) cut was removed from the loop. However, 
our power counting gives that it is exactly this cut that drives the enhancement of the one 
loop diagrams compared to the two-loop diagrams. 

Since the loop of type (f) gives the wave function at the origin in r-space, the formalism 
applied in Refs. [32] [33] [33] is basically identical to that used in the classic calculations for 
the decay of positronium into two photons. However, as discussed in detail in Ref. |36j . it is 
applicable only if the range of the transition potential from the constituents to the final state 
is significantly shorter ranged than the potential that formed the molecule — a scale not to 
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be mixed up with the size of the molecule which can be very large for a shallow bound state. 
However, the range of the binding momentum of the Z\> states is not known and might well be 
of the order of the range of the transition potential (at least as long as the final bottomonium 
is not a ground state). In such a situation in Ref. [36] it is proposed to calculate the full 
loop function for the transitions, as done here in our work, which in effect means to expand 
around the limit of a zero range potential that forms the bound state. In that paper it is 
also shown that the potentially most important corrections to the transition rate cancel, such 
that the uncertainty of the procedure is given by the binding momentum of the molecule in 
units of the range of forces and not of the order of the final momenta in units of the range of 
forces. This gives an additional justification for the approach we are using. We now discuss 
the formalisms of Refs. [321 [33| 134} [35] in some more detail. 

In Ref. [32] an effective field theory called X-EFT is used. It is valid for hadronic molecules 
with small binding energies so that the pion mass and the heavy meson hyperfine splitting 
are hard scales. The decays of the Z^ can be represented by a bubble with two vertices, one 

connecting the B^B* to the Z^ states and the other is a local operator for the B^ B*Tr(bb) 
coupling. The coefficient of the local operator depends on the pion energy, and is obtained 
by matching to the tree-level diagrams in heavy hadron chiral perturbation theory. For more 
details, we refer to Refs. [371 129]. The ratios of T(Z' T(3S»/T(Z T(3S» and 
T(Z' —7- hb(2P)ir)/T(Z — > hb(2P)ir) were calculated in Ref. [32] assuming z = z' . As outlined 
in the discussion below Eq. (|19p . in this case these ratios are just the ratios of phase spaces 
and thus our results for them agree for z = z' . The method of Ref. [33] is a phenomenological 
variant of the approach outlined above. 

Also in Ref. |34j the transition from the intermediate i?(*).B*-system to the final 7r(66) sys- 
tem was assumed to be local, however, here the strength of this local operator was calculated 
differently: the authors estimate it via the overlap integral of the bb component in the B^B* 
wave function with the outgoing bb pair in the presence of a dipole operator. While this 
procedure is certainly justified when there are IP states in the final state — here the relative 
momenta between the two B mesons are beyond 2.5 GeV and indeed in this case our effective 
field theory does not converge anymore (c.f. Sec. [2|) — we regard it as questionable for the 2P 
states. There is one more difference, namely the fact that in the formalism of Ref. [34] . the 
transitions to the JXbj(nP) final states are disconnected from those to the Trhb(nP) states, 
while in our approach they are connected as discussed in detail above. Thus, an experimental 
observation of the decay of one of the Z^ states to, say, 7Xtj(2P) would allow one to decide 
on the applicability of our approach. 

Similar to our work, in Ref. [35] the full heavy meson loop is evaluated, but regularized with 
a form factor. Absolute predictions are given for the transitions using a model to estimate the 
J3W B\*) (pb) coupling — it is difficult to judge the uncertainty induced by this. In a first step 
in that work a cut-off parameter was adjusted to reproduce each individual transition. It is 
found that the cut-off parameters needed for the hb~n transitions are typically larger and closer 
together than those needed for the Tir transitions. This hints at form factor effects being not 
very significant in the former decays. This interpretation is also supported by the observation 
that the ratios of decay rates — the same quantities as investigated here — are found to 
be basically independent of the form factor. In this sense the phenomenological studies of 
Ref. |35] provide additional support for the effective field theory calculation presented here, 
although a well-controlled error estimate cannot be expected from such a method. 
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5 Summary 



In this paper, we assume that both Zb and Z' h are hadronic molecules predominantly coupling 
to BB* and B*B*, respectively, in line with the data by the Belle Collaboration. As a 
consequence of this assumption, the Zb states can only couple through i?(*).B*-loops. Using 
NREFT power counting we argue that 

• the decay channels Z^' — > T(nS)ir as well as the transitions into the ground state P- 
wave bottomonia in the final state can not be controlled within the effective field theory, 
since higher loop contributions are expected to dominate the transitions; 

• model-independent predictions can be provided for Z^ — > hb{2P)ir and radiative decays 
Z^ -> X bj(2Ph- 

The ratios for Zb and Z' b decays into the same final states hb(mP)ir are consistent with the 
experimental data. Our results reflect the fact that those ratios are essentially the ratio of 
the corresponding phase space factors times the ratio of the Z^BB* and Z' b B*B* couplings 
squared. If further experimental analysis with higher statistics could underpin this fact, it 
would be a very strong evidence for the molecular interpretation since such a relation cannot 
be obtained from, e.g., a tetraquark structure. 

Furthermore, we calculate branching fractions for the final states X6j( m -P)7- They are 
predicted to be of order 10 -4 ~ 10~ 3 . Although this is clearly a challenge to experimentalists, 
a confirmation of these rates would strongly support the molecular picture. It is noted that 
the ratios of a certain Zb into bottomonia in the same spin multiplet are insensitive to the 
structure of the Zb, and may be obtained solely based on heavy quark spin symmetry. 
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A Nonrelativistic approach 

The basic three-point loop function worked out using dimensional regularization in D = 4 is 

/(rai,m 2 ,m 3 ,g) 

d D l I 1 1 



(2ir) D 

/^12^23_L 

16vr 



p_JL + ie 
mi 



tan 



+ tan 



l° + b 



12 



^23 



2a + c-c' 



2 v / ^(c r 



a — ie) 



(A.l) 



where nrii{i = 1,2,3) are the masses of the particles in the loop, /Xj,- = m^mj / '(rrii + rrij) are 
the reduced masses, 612 = m\ + 1712 — M, 623 = m 2 + ^3 + q° — M, with M the mass of the 
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initial particle, and 

2 



ft 



^23 \ -.2 oh ' o k i ^ 23 -2 

} q , c = 2^12612, c = 2^23»23 H q ■ 

m 3 J m 3 



In terms of the loop function given above, the amplitudes for Z£~ and Z, decays into hi,ir 



arc 



^+ft 6 = 2 ^ 9lZl sjM hb M Zh e m ^ z E k hh [I (M B , M B * , Ms* , q) + I(M B * , M B , M B * , 9)] 



(A.2) 



and 



= 2 ^ 9lZ2 JM hb M z ,e ijk q% l e k hb [I(M B . , M B * ,M B ,q) + I(M B * , M B * , M B * , 9)] , 

(A.3) 

respectively. In all these amplitudes, both the neutral and charged bottom and anti-bottom 
mesons have been taken into account. The amplitudes for Z^ into XfeJ7 read 

Mxmi = ~ \l\ i ^9iz^M Xb M Zb e tjk q i e 1 Zb e k 1 [I(M B , M B » , M B * , q) - 3 1{M B * , M B , M B , g)] , 

(A.4) 

A z?x*i = ~ 2i^(3egiz'v/M Xb M Zb e llk q l e j Zb e k I(M B * , M B * , M s * , q), (A.5) 
^o Xw7 = 2iPeg l z^/M Xb M Zb {q i gi k - q k g^)e Zb e^e k Xb I{M B , , Mb , M B * , g) (A.6) 
A z?w = " UPegiZy/M^MzMgi* - q k g^)e i Zb ei f e k b I(M B * , M B * , M B ,q), (A.7) 

and 

*4 z o Xb27 = V^iPeg^y/M^Mz^i^eiki + ^ ifem )4 6 44™/(M s * , M B * , M B * , g), (A.8) 
^o, XM7 = V2iPe 9l z^M Xb M Zb [q i (g jm e iM + <A^ m ) + q l e^ m + g"V^ 

x j z e k e l ™I(M B * , M B * , M B * , q), (A.9) 

respectively. 

B Relativist ic approach 

In this appendix, we formulate a Lorentz covariant framework for a parallel study of the 
Zj decays. In such a framework, more terms will appear as relativistic corrections that are 
generally neglected in the heavy quark limit. Therefore, a comparison between these two 
prescriptions will serve as a cross-check of the NREFT results, and more importantly as a 
confirmation of the validity of the NREFT power counting. 
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Hi 



The bottom meson fields are denned as 



2 



1 + ff 1 =7°i? 1 t 7° ; 



#2 = [P M V-75P] 



1 - f 



H =7 



P*V + 75^ 



7 °, F 2 = 7 °^ 7 °, 



where P* and P represent the (B* + , B*° , B*°) and (B + ,B°,B®) fields, respectively, which 
annihilate the corresponding particles, while P* and P are the fields of their antiparticles. 
The fields annihilating the S- and P-wave bottomonia are given by 



R 



dm 
r bb 



i + j> 

2 

1 + 1) 



(T^ 7/ , - %75 ) 



1 - f 



l-f 



respectively. The Lagrangians for S*-wave (P-wave) quarkonia and a pair of heavy mesons are 



C 



SBB 



C 



PBB 



ig 2 Tr 
%Tr 



R hl H 2al ^ M P la 



PjiH 2a ^H la 



+ H.c. , 
f H.c. 



Under the similar convention, the Z b field can be expressed as 



^75 



1 - j> 



(B.10) 
(B.ll) 

(B.12) 



and the effective interaction between the Z b and a pair of bottomed meson reads 



C 



ZBB 



iz (/) Tr 



The Lagrangian for the pion coupling to a pair of bottom mesons is [38, 39, 40J 

£ = igTt[H bW A^ a H a ] , 

where = (^d^ — £,d^)/2, with £ = ex.p(i\/2(f) / F w ) . The relativistic form of the La- 
grangian for the photon coupling to the bottom mesons is [30] 



+ H.c. , 



(B.13) 
(B.14) 



£ 7 = & I^ F ^ [H l^ vHa ] + ^F^[HlH a a, v }. 

With these effective Lagrangians, we obtain for the amplitudes of — > h^ir 
2y/2gg 1 z 



(B.15) 



Az h h h 



A Z!h 



2\/2gg 1 z 



e qV e hh e z ACo[M B ,M § *,M B *}+Co[M B *,M B ,M B *}) , (B.16) 



b"b 



(Co [M B * ,M B ,,M B *] + C [M B ,,M B *,M B ]). (B. 17) 



and to Y-7T 



2V2gg 



2Z 



(-er-e z \ q\ 2 (C [M B , M B * , M B *\ + 2C 2 [M B , M B * , M B * ] 



+ Co [M B * ,M B ,M B *} + 2C 2 [M b * , M B , M B * ] ) + g • e T q ■ e z (C [Mb* , M 5 , M B ] 
+ 2C 2 [M B *,M B ,M B )-C [M B *,M B ,M B *}-2C 2 [M B *,M B ,M B *))) , (B.18) 
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A Z , T = 2 ^ 92Z (e z ■ ey\q\ 2 (C [M B * , M g . , M B ] + C 2 [M B * , Mg* , M B ] 

+ C a [M B , , M B , , M B *} + C 2 [M B * , Mg* , M B *]) + q ■ e T q • e z {C [M B * , Mg, , M B ] 

+ C 2 [M B *,Mg*,M B ] - C [M B *,Mg*,M B *} - C 2 [M B * , Mg* , M B * ] ) ) . (B.19) 

The amplitudes for — > Xbol are 

A z o Xwl = i ^jf 1 e qvtieZb (3C Q [M B * , M B ,M B ] - C Q [M B ,Mg* , M B *]) (B.20) 

yl z o Xw7 = i2eg 1 zP(e Xbl -e J q-e z -q-e Xbl e 7 -e z )Co[M B *,Mg,M B *] (B.21) 

A z o X627 = i2v^e^^/3eSe a ^e XMffa C [M B ,M B -.,M B .] , (B.22) 

and 

A z?xwi = -^ieg lZ '(3e qve ^ z C [M B ,,Mg,,M B ,] (B.23) 

^^°X6i7 = ~ i2e 9i z ' P( e Xbi • e 7<? • - 9 • exw e 7 ' e z )C [M B *,Mg* } M B *] (B.24) 

A z ,o XM7 = i2v^e 5 i^(g a e ff „ £7£z + e°e^ £z )xS 2 aC'o[MB.,M B -.,M B .] , (B.25) 

respectively. Here, Cq is the standard relativistic three-point scalar function which is similar 
to the definition of / in App. |A] 



C [p 2 , {-qf, (p- g) 2 ,m|,mf,mg] 



(l2-m?)((p-Z)2-m!)((J-g)2-m§) 



2^, ' 



(B.26) 



where the incoming four-momentum of the Zf, is p, and the light outgoing particle four- 
momentum is q. For simplicity, we do not explicitly include the first three arguments in Co 
in the above formulaes. With the same convention, C 2 is a coefficient that arises from the 
tensor reduction of the three-point vector loop, 

c , s (2^) 4 - D f l"d D l 



{I 2 - m 2 )((p - If - m 2 )({l - qf - m 2 ) 
C lP » + C 2 {p-qY. (B.27) 



which is similar to the integral I± defined in Ref. [16J. 

To compare with the NREFT method, we note the following treatments in the heavy 
quark limit: 

• By expressing the heavy meson momentum Mqv^ = rriQV^ + kn in the heavy quark 
limit, where v is the heavy meson velocity and k is the residual momentum of the order 
of Aqcd, the derivative in Eq. (IB. 10}) only gives the difference of the residual momentum 
between the intermediate bottomed mesons [H]. Namely, the residual momentum k is 
the integral momentum in the meson loops instead of v. 

• The denominators (I 2 — m 2 )((p — I) 2 — m 2 )((l — q) 2 — m 2 ) of the meson loops contain 
two independent external momenta p and q. For the incoming momentum p, we set 
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p = Mzv following the convention of Ref. [16], where v = (1,0,0,0) defines the rest- 
frame of the initial particle. By doing this, the vector three-point function defined in 
the loop integrals becomes = Mzv^C\ + {M^v^ — q^)C2 which can be compared 
with the NREFT amplitudes term by term. It is interesting to note that the coefficient 
of C\ vanishes in the decay amplitude because of cancellation. 

• Another difference between this scheme and the NREFT comes from the contraction 
term q-er, where q is the external light meson (e.g. pion) momentum. Note that the time 
component appears as an additional contribution compared to the NREFT formalism. It 
is proportional to \q\E n /mY , which is a relativistic correction and relatively suppressed 
with respect to the space component \q\Ey/m^. 

• The relativistic corrections also arise from the mass difference between and tub*- 
One notices that in Eqs. (IB. 18ft and (IB, 19ft . the terms proportional to q ■ eyq ■ ez are 
given by the D-wave transition. In the heavy quark limit with ms = itib*, exact 
cancelations occur within the integral functions Cq and C2, respectively. This means if 
the Zh and Z' b are indeed the BB* and B*B* molecular states, respectively, their D- 
wave decays into the Ttt will be highly suppressed. This can be understood by noticing 
that the heavy quark spin decouples from the system in the heavy quark limit, and 
the total angular momentum of the light quarks s qq - is a good quantum number. If the 

Z^p states are S'-wave B^B* hadronic molecules as assumed here, there is no spatial 
angular momentum in the system. Thus, the light quark system has s qq - = or 1, and 
therefore cannot couple to a spinless pion in a D-wave. A similar statement was made 
very recently in Ref. [52] . Here we notice that the decay of a bbqq tetraquark state with 
= 2~ would decay into bbn dominantly in a D-wave. 

• We also note the convention for the sums of the polarizations for the vector and tensor 
particles: 4(^)4^) = + = and e^At = +9^9 ua ) ~ 
^gVg^. They allow us to separate out the relativistic contributions in the scalar loop 
diagrams which can then be compared with the NREFT formulaes explicitly. 

The relativistic formalism is advantageous for providing a cross check of the NREFT 
results and singling out the effects arising from the relativistic corrections. For the bottomo- 
nium and bottomed meson system discussed here, it shows that the relativistic corrections are 
indeed small and these two methods are consistent with each other. Alternatively, it should 
be cautioned that in the relativistic formalism, when the terms proportional to {qlE^/rnqq 
are not obviously suppressed in comparison with the three-momentum \q\, i.e. E n /m,QQ is 
sizeable, a different scheme with a form factor [43 \ HJJ might be used to control the large 
momentum contributions and non-local effects of each couplings. 
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